Layer-thickness effects in quai-two dimensional electron liquids. 



in 
o 
o 

(N 



M.W.C. Dharma-wardan^ 

Institute of Microstructural Sciences, 
National Research Council of Canada, 
Ottawa, Canada. KIA 0R6 

(Dated: February 2, 2008) 

We use a mapping of the quasi-2D electron liquid to a classical fluid and use the hypernetted-chain 
equation inclusive of bridge corrections, i.e., CHNC, to calculate the electron pair-distribution func- 
tions and exchange-correlation energies as a function of the density, layer width, spin-polarization 
and temperature. The theory is free of adjustable parameters and is in good accord with recent 
effective-mass and spin-susceptibility results for HIGFET 2-D layers. 
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The 2-D electron systems (2DES) present in GaAs 
or Si/Si02 nanostructures access a wide range of elec- 
tron densities under controlled conditions, providing a 
wealth of information 1] which is of basic and technolog- 
ical importance. The 2DES is in the x-y plane and also 
has a transverse extension in the lowest sub-band of the 
hetero-structure0- The higher subbands are far above 
the Fermi energy and no upward excitations are possible. 
Although the z-motion is confined to the lowest subband, 
realistic layers may have widths of ~ 600 A or more, and 
this corresponds to ~ 6 effective atomic units in GaAs. 
Recent experiments and theory have focused on these 
layer-thickness effects 0, H [H EH- The physics of 
such quasi-2DES depends on the density parameter r^, 
the layer thickness w which labels the z-charge distribu- 
tion, the spin-polarization C, and the temperature T. The 
2D density n defines the mean-disk radius ~ (vrn)^^/^ 
per electron, expressed in effective atomic units which 
depend on the bandstructure mass and the "back- 
ground" dielectric constant £(,. Although is the "small 
parameter" in perturbation theory (PT), here it is simply 
the electron-disk radius and PT is not used. 

Finite-width effects of the 2DES arise also in the quan- 
tum Hall effecting, where an "unperturbcd-g" approxi- 
mation, which uses the pair-distribution function (PDF) 
of the ideal 2DES and the quasi-Coulomb potential W{r) 
of the thick 2D layer are used to calculate energies. 

While diagrammatic methods propose "turn-key" pro- 
cedures for calculating many-body properties, they work 
only for weakly coupled (small r^) systems, where RPA- 
like approximations may be used. Varied results can 
be obtained using various approximations which go "be- 
yond" RPA. Unfortunately, an approximation success- 
ful with one property usually fails for other properties. 
Such methods lead to negative PDFs, incorrect local- 
field corrections in the response functions, and disagree- 
ment with the compressibility sum rule etc., and in- 
correct predictions of spin-phase transitions (SPT) at 
quite high densities. In fact, alternative approaches were 
sought by Singwi, Tosi et al. (STLS)'T^, and Ichimaru 



et al.[l2j, and also within the Feenberg-type correlated- 
wavefunction methods [ij]. Most of the currently avail- 
able results for strongly-coupled systems have been 
generated using correlated-wavefunction approaches via 
Quantum Monte Carlo (QMC) simulations. However, 
QMC remains a strongly computer-intensive numerical 
method which is best suited for the study of simple 
"bench-mark" systems. 

We have recently introduced a conceptually and 
numerically simple, adequately accurate method for 
strongly-correlated quantum systems at zero and finite 
temperatures, where the objective is to work with the 
PDF of the quantum fluid, generated from a classical 
Coulomb fluid whose temperature Tq is chosen to re- 
produce the correlation energy of the original quantum 
fluid at T = 0. The classical PDFs are obtained from 
an integral equation equivalent to a classical Kohn-Sham 
equation where the correlation effects are captured as a 
sum of hyper-netted-chain (HNC) diagrams and bridge 
diagrams. This method, known as a Classical mapping 
to an HNC form, i.e, CHNC, was applied to the 3D 
and 2D electron fluids[ll [H [H [13, to dense hydro- 
gen fluid ligll. and also to the two-valley system in Si- 
M0SFETs|2a E^ . In each case we showed that the 
PDFs, energies, etc., obtained from CHNC were in ex- 
cellent agreement with comparable QMC results [l^, In- 
variants of the method have also been discussed by Bu- 
lutay and Tanatar, and by Khanh and Totsuji(25j]. 

The advantage of CHNC is that it affords a sim- 
ple, semi-analytic theory for strongly correlated systems 
where QMC becomes prohibitive or technically impos- 
sible to carry out. The classical-fluid model allows for 
physically motivated treatments of complex issues like 
three-body clustering etc., via the statistical mechan- 
ics of hard-disk reference fiuidsjl^. The disadvantage 
of the method, typical of such many-body approaches, 
is that it is currently an "extrapolation" method tak- 
ing off from the results of a model fluid. For 2D sys- 
tems, the fully spin-polarized ideally-thin uniform fluid 
is the model fluid[r3]- As the CHNC method has been 
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described in previous work, we do not give a detailed 
account here. This study is a simple, direct application 
of CHNC to the quasi-2D potential W{r) = V{r)F{r), 
where V{r) = 1/r and F(r) is a form-factor accounting 
for the modifications arising from the thickness effect. 

Most of our calculations are for the Fang-Howard(FH) 
approximation 2J to the charge density n{z) contained in 
HIGFET structures of the type used by Zhu et al|l|- In 
this case, F{k), the form factor in fc-space, has an ana- 
lytic form but F(r) is numerically determined. It is tech- 
nically convenient to work with an equivalent constant- 
density model (CDM) for which analytic forms are avail- 
able for F{r) as well as F{k). We present a potential 
W{r, w) for a CDM which is electrostatically equivalent to 
the the 2D potential for any n{z), e.g., the FH-potential 
defined by the parameter b. The method of replacing an 
inhomogeneous distribution by a uniform distribution is 
suggested by the observation that the non-interacting to- 
tal correlation function h'^{r) = g^{r) — 1 has the form 
~ nirY' , where n[r) is the density-profile around the 
Fermi hole. We replace the inhomogeneous n{z) by a 
slab of constant-density ricd which is easily shown to ac- 
curately recover the electrostatic potential of n{z) in the 
2-D plane. 

Tied = 1/w = j n{z)'^dz (1) 

Since the subband distribution is normalized to unity, 
the width w of the CDM is simply l/ricd- Hence quasi- 
2D layers can be labeled by their effective width w. The 
CDM width w for the n{z) labeled by the Fang-Howard 
h is w = 16/(36), and differs from the commonly quoted 
"thickness", 3/6. This provides an explicit example of 
the replacement of an inhomogeneous distribution by a 
homogeneous distribution, already suggested in Ref. j27l | 
and used for 2-clcctron atoms. The quasi-2D potential 
for a CDM of width w is given by 

W{r) = V{r)F{s), s^r/w, t^y/{l + s^) 
F{r) = 2s[logit /s)+t] (2) 

This tends to 1/r for large r. The short-range behaviour 
is logarithmic, and weaker than the 1/r potential. The 
form factor F{r) for a HIGFET with r^ = 5 is shown 
in the inset to Fig. ^ The /c-space form of the CDM 
potential is: 

W{k,w) = V{k)F{p), p^kw (3) 
F{p) = i2/p){{e-P^l)/p+l} (4) 

The form factors F{s) and F{p) tend to unity as w ^ 0. 
These r-space and fc-space analytic forms of the CDM 
indirectly lead to analytic formulae for the FH from. 
The W{r) of the CDM is equivalent to that from the 
original inhomogeneous distribution, and only W{r) en- 
ters into the exchange-correlation and g{r) calculations. 
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FIG. 1: (a)The effective mass m* of an ideal 2D layer (w = 
0) obtained here are compared with the the QMC data of 
Ref. I31II and the G+&lG- / D calculation of Asgari et al[lo(|. 
(b)The effective mass uih in a HIGFET using the CHNC 
ideal-2D finite-T g{r) in the "unperturbed-g" approximation, 
i.e., Eq. 121 The HIGFET calculation of Asgari et al. is also 
shown. 



In the case of GaAs-HIGFETS, if the depletion density 

could be neglected ^9], r^ specifies the h parameter and 

hence the width w of the CDM. Then 6^ = 33/(2r2) and 
2/3 

w = 2.09494rs . When the exchange-correlation energy 
Exc{b,rs) is included, 6 changes to 6*. The correction is 
~ 2 — 3% at low r^ ^ 1 and decreases as r^ increases. We 
have used b* in our fully self-consistent CHNC calcula- 
tions. 

Exchange and correlation in quasi-2D layers. — The 
noninteracting-2D correlation function, h^{r) — g''\r) — l, 
yields exact exchange energies for arbitrary w as well 
as for w — Q (the ideal 2DES). On the other hand, 
the correlation-energy evaluation needs the PDFs of the 
quasi-2D potential AVF(r, w), at many values of the cou- 
pling constant A. These g{r,(^,w,X) can be calculated 
using the CHNC. However, the use of the unperturbed-g 
approximation used in Quantum Hall studies can be 
useful here too. De Palo et al.^ have in fact exploited 
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such an approach where the g{r, (,w = 0) of the ideally 
thin layer are used to calculate a correction energy A 
given by, 

A = (n/2) J 2nrdr[W{r) ~ V{r)]h{r, C,w^Q) (5) 

Then the total ExcirsXiw) is obtained by adding to A 
the known Exc of the ideally-thin system. The "unper- 
turbed" g{r) needed in Eq. [H at T = are the QMC 
5(r)^. Equation El can also be applied to the finite- 
T ideal g{rs, Ci T) obtainable from CHNC. However, this 
approach is found to be insufficient for calculating m* . 

De Palo et al0 have performed Diffusion Monte Carlo 
simulations at =5 for HIGFETS with b = 0.8707, i.e, 
a CDM width i(j=6.1256 a.u., and find that the error 
compared to the full simulation is about 2%. Since the 
ratio w/ts decreases in the HIGFET as increases, the 
HIGFET approaches the thin-layer model for large Tg. 
Hence Eq. [Sjis satisfactory for > 5, and unreliable for 
small Ts, e.g., below = 3. Equation neglects the 
renormalization of the kinetic energy, correction of b to 
6*, as well as the changes in g{r) due to the changed 
potential. We have used both the full CHNC which ac- 
counts for all these effects, and also the "unperturbed-^" 
approximation, Eq. |S1 and find that the latter is indeed 
satisfactory for the calculation of Exc and the suscepti- 
bility enhancement m*g* = Xs/xpi where Xs, Xp are the 
interacting and Pauli spin susceptibilities. 

Correlation energy at finite temperatures. — The corre- 
lation contribution to the Helmholtz free energy of ideal 
{w — 0) or thick layers [w > 0) is readily calculated via 
Eq. using the CHNC-finite-T g{r). For example, con- 
sider r^ = 5, and T/Ef=0-05, with 77 = 0.3718, i.e., the 
packing fraction of the hard-disk fluid used to mimic the 
three-body terms in the extended HNG equation. Then 
Exc for the ideal 2D at C = is -0.16902 a.u., while the 
HIGFET Exc from Eq. and from the full CHNC are - 
0.11197 a.u. and -0.11467 a.u. respectively. As discussed 
in earlier work |2lll2^ . the exchange free energy and 
the correlation free energy Fc at very low T contain log- 
arithmic terms which cancel with each other, so that the 
sum Fxc = Fx + Fc is free of such terms. At = 5 
the cancellation is good to about 75%, and this improves 
as rs increases. Although the two-component fluid (up 
and down spins) involves three distribution functions, 
we have, as before ^3 1 used only one hard-disk bridge 
function, B12, as clustering effects in gu are mostly sup- 
pressed by the Pauli-exclusion. This is a satisfactory ap- 
proximation for Ts > 3, and for smaller rs B12 becomes 
comparable to Bu, but negligible. Also, B12 is assumed 
to be independent of C and identical to the ideal- 2D 
{w = 0) bridge function. 

The spin- susceptibility, m* and g* . — The results for 
the Fxc{rs,C,T,w) of the ideal or thick 2DES suffice to 
calculate the spin-susceptibility enhancement, the effec- 
tive mass m* and the effective Lande factor g* . We cal- 
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FIG. 2: (a)The spin-averaged g(r) at = 5 and T/Ef=0.1, 
for an ideally thin layer, and for a HIGFET calculated us- 
ing CHNC. (b)The effective mass in a HIGFET using 
the fully self-consistent CHNC calculations for the energy at 
finite-T. The experimental data are from Tan et alQ- 



culate the quantity A= 1 + B{x), where B{x) is the ra- 
tio of the second derivative with respect to x of Fxc{x), 
and ^0(2;), where x is C for the spin susceptibility cal- 
culation, and X = T for the m* calculation. Fq[x) is 
the noninteracting Helmholtz free energy. The CHNC 
is used to obtain any data (e.g, at finite -T) unavail- 
able from QMG. The effective mass obtained for the 
HIGFET using Eq. where the ideal 2D-finite-r g{r) 
is used, is presented in the lower panel of Fig. ^ The 
upper panel shows the ideal 2D-layer m* , in good agree- 
ment with the four QMC values. This contrasts the m* 
proposed by Asgari et al, denoted G+hG- / D in their 
paper A study of the local-field factor of the 2D re- 
sponse function[2J| shows that the formation of singlet- 
pair correlations is complete by ~ 5, and after that the 
structure of the fluid remains more or less unchanged, un- 
til the spin-phase transition (SPT) is reached. The rapid 
rise in to* before ^ 5 and the subsequent slow-down 
are probably related to the formation and persistence of 
singlet structure in the 2D fluid. 
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FIG. 3: The Lande g-factor for the ideal 2D Layer is ob- 
tained from the QMC rn g* of ref. 33], divided by the CHNC 
ideal 2D m*. The experimental HIGFET g\i (boxes) is from 
the m* g* of Zhu et al., divided by the m* of Tan et al. We 
also show g*ii calculated from the m* g* and the mjf results 
for the HIGFET. The inset shows the spin-susceptibility en- 
hancement m* g* from the C, dependent energies calculated 
from Eq.|3 where the ideal 2D g(r) is used, and from the full 
CHNC calculation using the gir) consistent with the quasi-2D 
potential. 



In the lower panel of Fig. ^ we present the effec- 
tive mass of the electrons in the HIGFET. The 
"unperturbed-^" model, Eq.jS] is unable to approach the 
experimental results (Fig. of Tan et al. Hence we re- 
peat the m* calculations using the fully self-consistent 
g{r) from CHNC. The g(r) for the ideal 2DES and 
the fully self-consistent g(r) for a HIGFET at — 5, 
T/Ef — 0.1 are shown in the upper panel of Fig. |21 
The weakened repulsion in W{r) for small-r in the quasi- 
2DES allows easier electron clustering, strongly boosting 
the g{r) near r = 0. The HIGFET m*, evaluated from 
the second derivative of Fxc{T) in the range T/Ep — 0.05 
to 0.1, is shown in the lower panel. Note that m* in our 
work increases from the ideal [w = 0) system to the 
HIGFET, where as Asgari et al. predict the opposite. 

Enhancement of the spin susceptibility. De Palo et 
have calculated m*g* from the ideal QMC g{r), us- 
ing the "unperturbed-^" approximation, and shown that 
they obtain close ag reement with the data for very nar- 
row 2-D systems |33| as well as for the thicker 2DES in 



HIGFETSi, as shown in the inset to Fig. |3 We ob- 
tain results in close agreement when the "unperturbed- 
5" approximation, Eq. 5, is used. The m*g* from the 
fully self-consistent g{r) are also shown in the inset, and 
shows a deviation of ~ 5%. This deviation is probably 
a short-coming of the CHNC model, associated with the 
use of a spin-independent bridge function. 

The effective Landee-g factor. — The agreement be- 
tween the experimental Xs/xp ratio (i.e, m*g*) and 
the theoretical results, especially those from QMC for 
the ideal 2DES and the HIGFET suggests that m*g* 
is known with some confidence. Hence we may extract 
the effective Lande-5 factors for the ideal 2DES and the 
HIGFET (see Fig. ISJ, using the available m* values. 
The strong increase in m* with Ts in HIGFETS implies 
that the g* -factor is less sensitive to r^. This is simi- 
lar to the behaviour of the two-valley system found in 
Si-MOSFETS where there is no SPT^ The ad- 

ditional inter-valley Coulomb interactions in MOSFETS 
weakens the role of exchange and enhanced the effect of 
singlet-pair cluster effects, leading to a strong increase 
around = 5. In the HIGFET, the the weakening of 
the repulsive interaction increases clustering, and boosts 
m*, since quasi-particles no longer move freely, but have 
to drag a cloud of electrons associated with the enhanced 
g{r) near r/rg 0. This shows that the m* calculation is 
very sensitive to the accurate evaluation of the "on-top" 
value of g{r = 0). In the CHNC, clustering and such ef- 
fects are controlled by the hard-disk bridge function, and 
the diffraction correction (de Brogiie momentum kth JJi] ) 
used to describe the quantum-scattering of two electrons. 
In this work we have simply used the 77 and kth of the ideal 
2DES, i.e., w — 0. Better agreement of m* with experi- 
ment would require an evaluation of these as a function 
of w. Similarly, complete QMC runs for finite w would 
require back-flow and three-body functions different from 
the usual RPA-like treatment of the 2DES. 

Conclusion- We have used a single theoretical frame- 
work, i.e., the CHNC, with no parameters other than 
those previously used for the ideal 2D svstem[l7||. to cal- 
culate the Fxc{rsTC,,T^w)^ and hence the m* and g* of 
thick 2D layers. The results suggest that the enhanced 
mass in HIGFETS arises from the strong short-ranged 
correlations created by the weakened Coulomb repulsion 
due to the thickness effect. 
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